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Abstract 

Background: There is increasing appreciation for sexually dimorphic effects, but the molecular mechanisms 
underlying these effects are only partially understood. In the present study, we explored transcriptomics and 
epigenetic differences in the small intestine and colon of prepubescent male and female mice. In addition, the 
microbiota composition of the colonic luminal content has been examined. 

Methods: At postnatal day 14, male and female C57BL/6 mice were sacrificed and the small intestine, colon and 
content of luminal colon were isolated. Gene expression of both segments of the intestine was analysed by 
microarray analysis. DNA methylation of the promoter regions of selected sexually dimorphic genes was examined 
by pyrosequencing. Composition of the microbiota was explored by deep sequencing. 

Results: Sexually dimorphic genes were observed in both segments of the intestine of 2-week-old mouse pups, with a 
stronger effect in the small intestine. Amongst the total of 349 genes displaying a sexually dimorphic effect in the small 
intestine and/or colon, several candidates exhibited a previously established function in the intestine (i.e. Nts, Nucb2, 
Alox5ap and Retnly). In addition, differential expression of genes linked to intestinal bowel disease (i.e. Ccr3, Cell 1 and Tnfr) 
and colorectal cancer development (i.e. Wtl and Mmp25) was observed between males and females. Amongst the genes 
displaying significant sexually dimorphic expression, nine genes were histone-modifying enzymes, suggesting that 
epigenetic mechanisms might be a potential underlying regulatory mechanism. However, our results reveal no significant 
changes in DNA methylation of analysed CpGs within the selected differentially expressed genes. With respect to the 
bacterial community composition in the colon, a dominant effect of litter origin was found but no significant sex effect 
was detected. However, a sex effect on the dominance of specific taxa was observed. 

Conclusions: This study reveals molecular dissimilarities between males and females in the small intestine and colon of 
prepubescent mice, which might underlie differences in physiological functioning and in disease predisposition in the 
two sexes. 

Keywords: Small intestine, Colon, Sexually dimorphism, Gene expression, Microbiota, Histone modification, DNA 
methylation, Epigenetics, Chromosomes, Prepubescent 
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Background 

There is growing recognition that males and females 
differ with respect to basic physiology, body compos- 
ition and susceptibility to and progression of a broad 
variety of non-communicable diseases, as well as in the 
response to pharmacological treatment. The molecular 
mechanisms underlying these sexually dimorphic effects 
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are currently largely unresolved. Increased knowledge on 
these mechanisms might contribute significantly to dis- 
ease prevention and treatment, for instance by optimizing 
dietary recommendations and pharmacological protocols 
in a sex-specific way. 

Differences between males and females not only become 
clearly apparent post-puberty in response to increasing 
circulating hormone levels but are also present at a much 
younger age. During the initial phase of embryonic devel- 
opment, differences between males and females are 
already detectable [1-4]. Immediately thereafter, male 
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expression of the Y chromosome Sry gene initiates 
molecular and cellular cascades in the undifferentiated 
gonads of males leading to differentiation of testes and 
divergent conversion between male and female devel- 
opment [5]. Insight into how the sex chromosomes 
contribute to the differences between males and females 
has strongly increased in the last decade. The male Y 
chromosome is relatively small (60 Mb) and contains only 
a limited number of genes (-60 in humans and 12 in 
mice) [6,7]. Even though these genes are exclusively 
expressed in males, a female variant for some of them has 
been identified [8]. The X chromosome, which is present 
as one copy in males and as two copies in females, is 
much larger in size (155 Mb) and contains a much larger 
number of genes (~ 1,500). To ensure dosage equivalence 
between males and females, one of the two X chromo- 
somes is transcriptionally silenced in females [9,10]. X in- 
activation silences randomly either the paternal or the 
maternal X chromosome, and Xist, a 17-kb spliced and 
polyadenylated RNA with no coding capacity, plays a 
crucial role in this process [11-14]. Some genes located on 
the X chromosome escape silencing in female somatic 
cells, including a number of genes containing a homologue 
on the Y chromosome such as the well-known escapee 
KdmSc [15,16]. 

Genome-wide gene expression analysis has revealed 
that sex differences occur in a wide variety of organs 
and tissues in humans as well as in rodents, including 
the liver, heart, kidney, brain, muscle, placenta, gonads 
and adipose tissue [17-25]. Sexually dimorphic genes have 
been best characterized in the liver [19,22-24,26-32]. Most 
studies analysing sexually dimorphic effects have been 
carried out in adult tissues where sex hormones likely 
contribute strongly to the observed effects. Conforto and 
Waxman have recently shown that— although by far the 
most pronounced difference in gene expression was found 
in the livers of adult mice— a significant sexually dimorphic 
expression was already detectable in the prepubescent 
mice [26]. Similar results were reported by Kwekel and col- 
leagues for the kidney in rats [21]. 

Until now, sexually dimorphic effects in the intestine 
have only marginally been described. Oestrogen-mediated 
health effects in the intestine have been reported [33-35], 
and sexually dimorphic gene expression in the small intes- 
tine (SI) has been shown in adult mice [36]. Information 
regarding sexually dimorphic expression in the SI at a pre- 
pubescent stage is still lacking, and no information is cur- 
rently available regarding sex-regulated differences in gene 
expression in the colon. While in the SI food is degraded 
by enzymatic digestion prior to absorption, the dominant 
process in the colon is fermentation by microorganisms. 
Inoculation and colonization of the newborn sterile colon 
starts immediately during and after delivery [37]. The ex- 
istence of a strong host-microbiota interaction is currently 



commonly acknowledged. Furthermore, the importance of 
the gut microbiota for the development of the metabolic 
and immune system early in life is nowadays commonly 
appreciated [37]. Since gene expression of the colon is 
directly linked to the microbiota composition of the colon 
lumen [38-40], it can be speculated that sexually 
dimorphic expression of genes in the colon affects the 
microbiota composition. 

Sex-related differences in gene transcription might be 
attributable to various mechanisms. It could relate to the 
location of genes on the sex chromosomes, as well as to 
epigenetic marks associated with the genes. Sex-specific 
DNA hypersensitive sites [41] and differences in DNA 
methylation [23] and in histone marks [22,42] have been 
reported between males and females. Whether these 
mechanisms function equally in different organs and are 
globally regulated or mediated at a gene-specific level 
needs to be further investigated. 

Tissue-specific sexually dimorphic gene expression is 
not well-characterized in either the SI or the colon, par- 
ticularly in prepubescent subjects. To address this gap, 
we performed differential gene expression profiling in 
the SI and colon of 2-week-old male and female mouse 
pups by applying microarray (MA) analysis. Pyrosequencing 
was performed for the genes exhibiting the strongest 
differential gene expression to measure whether DNA 
methylation differences are present in the promoter re- 
gions of the sexually dimorphic genes, thereby possibly 
affecting the expression state. The microbiota compos- 
ition of the colon luminal content was explored by deep 
sequencing and investigated for possible sex differences. 

Methods 

Ethics statement 

The national and institutional guidelines for the care 
and use of animals were followed, and the experimental 
procedures were reviewed and approved by the Ethics 
Committees for Animal Experiments of the University of 
Groningen, The Netherlands (Ethics registration code: 
5709). 

Animals and diets 

Female C57BL/6 mice (5 weeks of age) were purchased 
from Harlan (Horst, The Netherlands) and housed indi- 
vidually in the light- and temperature-controlled facility 
of the University Medical Center Groningen (lights on 
7:00 am-7:00 pm, 21°C). The mice had free access to 
drinking water and were fed a semi-synthetic low-fat 
control diet (3.85 kcal/g; 10 E% fat, 20 E% protein, 70 E% 
carbohydrate; D12450B, Research Diets, New Brunswick, 
NJ, USA). After 6 weeks (pre-treatment period), the female 
mice were mated with males and, in case conceiving failed, 
allowed to re-mate. Mice were allowed to deliver spontan- 
eously and were left undisturbed with their litters for 24 h. 
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Litter sizes were standardized to five to seven pups, to 
ensure no litter was nutritionally biased due to lower or 
higher litter size. By natural circumstances, the litter size of 
some dams was reduced further, but not changing the 
overall male/female ratio significantly. Throughout preg- 
nancy and lactation, the dams received the same diet. After 
2 weeks of lactation, the offspring were sacrificed by heart 
puncture under isoflurane anaesthesia. The SI and colon 
were isolated from each mouse, snap-frozen in liquid 
nitrogen and stored at -80°C until further use. In total, 
nine female and six male mice were included derived 
from three different litters. Physiological and molecular 
effects observed in the liver of these mice have been re- 
ported previously [43]. 

RNA isolation 

Total RNA was isolated from colon and SI samples as 
described previously [43]. In brief, TRIzol reagent 
(Invitrogen, Breda, The Netherlands) was used and the 
samples were treated with DNAse and purified on columns 
(RNeasy Micro Kit, Qiagen, Venlo, The Netherlands), 
all according to the manufacturers instructions. Puri- 
fied RNA was immediately stored at -80°C until further 
use. RNA concentrations were determined using the 
NanoDrop ND-1000 UV-vis spectrophotometer (Isogen, 
Maarsen, The Netherlands). RNA integrity was verified 
on an Agilent 2100 Bioanalyzer with the 6000 Nano Kit 
using the Eukaryote Total RNA Nano Assay according 
to the manufacturers instructions (Agilent Technolo- 
gies, Amsterdam, The Netherlands). Samples were con- 
sidered suitable for hybridization when they showed 
intact bands of 18S and 28S ribosomal RNA subunits, 
displayed no chromosomal peaks or RNA degradation 
products and had a RNA integrity number (RIN) 
above 8.0. 

Microarray hybridization and analysis 

Per offspring, colon and SI samples were analysed as 
described previously [43], and in total, six male and six 
female samples were used in this analysis (three of the 
nine females were not included in the microarray ana- 
lysis due to budget limitations). In brief, 100 ng of puri- 
fied RNA was used for the preparation of labelled cDNA, 
applying the Ambion Whole Transcript (WT) Expression 
Kit (Life Technologies, Carlsbad, CA, USA) in combin- 
ation with the Affymetrix GeneChip WT Terminal Label- 
ing Kit (Affymetrix, Santa Clara, CA, USA). All samples 
were hybridized at one time point to Affymetrix GeneChip 
Mouse Gene 1.1 ST arrays according to standard Affyme- 
trix protocols. Quality control and normalization were 
performed using Bioconductor software packages inte- 
grated in an on-line pipeline [44]. Normalized expression 
estimates of probe sets were computed by the robust 
multiarray (RMA) analysis algorithm available in the 



Bioconductor library AffyPLM using default settings [45] . 
Probe sets were redefined according to Dai et al. [46] and 
assigned to unique gene identifiers (IDs) of the Entrez 
Gene database, resulting in 21,187 assigned Entrez IDs. 
Array data were submitted to the Gene Expression Omni- 
bus and are available under accession number GSE57516. 

Bioinformatic analysis 

Of the 21,187 defined genes covered by the MA, only 
genes with an intensity value of >20 on at least five ar- 
rays, represented by at least seven probes per gene on the 
array and an interquartile range (IQR) >0.1 were selected 
for further analysis whereby the colon and SI were analysed 
separately. Pancreas-specific genes (http://biogps.org) were 
removed from the analysis in order to omit an effect of po- 
tential pancreatic contamination. For the SI in total 12,297 
and for the colon 14,400 genes were included in the subse- 
quent analysis; 10,958 genes met the above-described cri- 
teria in both the SI and colon. The top 1,000 most variable 
genes were used for principal component analysis (PCA) 
using MultiExperiment Viewer version 4.8.4 [47,48] with 
an eigenvalue of 1.0 as a cutoff for identification of contrib- 
uting components. Signal 2log ratios, which represent fold 
changes (FC), and related significances of change were cal- 
culated from the mean signal intensities of the female and 
male groups for analysis of the sex-specific differential gene 
expression using intensity-based moderated t statistics 
(IBMT) implementing empirical Bayes correction [49]. 
Resulting 2log ratios and p values were applied for further 
descriptive bioinformatic analysis of the data. Ingenuity 
pathway analysis (IPA; Ingenuity® Systems, www. ingenuity, 
com) was used to relate the MA data to networks of dis- 
eases and biological functions. Comparison of male and fe- 
male expression patterns in the SI and colon on genes 
located on the autosomes was carried out by generating heat 
maps using MultiExperiment Viewer version 4.8.1 [47,48]. 

DNA isolation from the SI and colon 

Genomic DNA was isolated from the SI and colon by 
using the DNeasy® Blood and Tissue Kit (Qiagen, Venlo, 
The Netherlands) as described before [50] according to the 
manufacturers instructions. The DNA was treated with 
RNase and eluted in Qiagen elution buffer AE. DNA purity 
and quantity were checked spectrophotometrically (ND- 
1000, NanoDrop Technologies, Wilmington, DE, USA). 

Bisulfite conversion and DNA methylation analysis 

For each sample, 500 ng of genomic DNA was bisulfite- 
treated using the EZ-96 DNA Methylation- Gold™ Kit 
(Zymo Research, Irvine, CA, USA) and eluted in 14 \A of 
M-Elution Buffer. DNA methylation analysis was per- 
formed using PyroMark™ pyrosequencing technology 
(Biotage AB, Uppsala, Sweden). Primers were designed 
using PyroMark software, and the sequences of the primers 
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used are listed in Additional file 1. The PCR reactions were 
performed in a total volume between 25 and 50 ul, and the 
volume of bisulfite-treated genomic DNA used was always 
1/20 of the total PCR volume. PyroMark PCR Master Mix 
and CoralLoad Concentrate were used according to the 
manufacturers instructions, and 0.2 \iM of each primer 
(Qiagen, Venlo, The Netherlands) was used. The following 
thermal cycling conditions were applied: 15 min at 95°C, 
followed by 45 cycles of 94°C for 30 s, tempX (gene- 
specific, see Additional file 1) for 30 s, and 72°C for 
40 s, followed by a final elongation step at 72°C for 
10 min. The PCR product was bound to Streptavidin 
Sepharose HP beads (GE Healthcare, Uppsala, Sweden) 
and purified and made single-stranded using the Pyro- 
sequencing Vacuum Prep Tool according to the manufac- 
turer s instructions (Qiagen, Venlo, The Netherlands). 
Sequencing primers (for sequences, see Additional file 1) 
were annealed to the purified single-stranded PCR prod- 
uct, and pyrosequencing was performed using the Q24 
Pyrosequencing System (Qiagen, Venlo, The Netherlands). 
CpG methylation was analysed with the provided software. 

Bacterial DNA extraction 

DNA from all six male and nine female sacrificed mice 
was extracted from the freeze-dried luminal content of 
the colon using the method described by Salonen et al. 
[51]. In short, approximately 0.1 g was used for mechan- 
ical and chemical lysis using 0.5 ml buffer (500 mM 
NaCl, 50 mM Tris-HCl (pH 8), 50 mM EDTA, 4% SDS) 
and 0.25 g of 0.1-mm zirconia beads and 3-mm glass 
beads. Nucleic acids were precipitated by addition of 
130 \A of 10 M ammonium acetate, using one volume of 
isopropanol. Subsequently, DNA pellets were washed 
with 70% ethanol. Further purification of DNA was per- 
formed using the QIAamp DNA Stool Mini Kit (Qiagen, 
Hilden, Germany). Finally, DNA was dissolved in 200 \A 
Tris/EDTA buffer, and its purity and quantity were 
checked spectrophotometrically (ND-1000, NanoDrop 
Technologies, Wilmington, DE, USA). 

Library preparation for 16S rRNA pyrosequencing 

Universal primers were applied (forward primer, 5'-CC 
ATCTCATCCCTGCGTGTCTCCGACTCAGNNNNNNA 
CTCCTACGGGAGGCAGCAG 3 ; reverse primer, 5 
CCTATCCCCTGTGTGCCTTGGCAGTCTCAGCRRCA 
CGAGCTGACGAC-3 ') for amplification of the V3-V6 
region of the 16S rRNA gene as described before [52]. 
The forward primer includes a sample-specific six-base 
barcode (NNNNNM) to tag each PCR product. The amp- 
lification PCR consists of 2 \A microbial genomic DNA, 
16 ul master mix (1 ul KOD Hot Start DNA Polymerase 
(1 U/ul; Novagen, Madison, WI, USA), 5 [A KOD buffer 
(lOx), 3 [A MgS0 4 (25 mM), 5 [A dNTP mix (2 mM 
each), 1 ul (10 uM) of each forward and reverse primer) 



and 32 \A sterile water (total volume 50 \A). PCR condi- 
tions were as follows: 95°C for 2 min followed by 35 cy- 
cles of 95°C for 20 s, 55°C for 10 s and 70°C for 15 s, 
ending with a last step of 72°C for 10 min to ensure 
complete amplification of the target region. From each 
sample, 5 ul was electrophoresed on a 1% agarose gel 
and the approximately 750-bp PCR amplicon was sub- 
sequently purified using the MSB Spin PCRapace kit 
(Invitek, Westburg, The Netherlands) followed by a sec- 
ond purification step using PureLink columns (Invitrogen, 
Breda, The Netherlands). The concentration was checked 
with a NanoDrop 2000 spectrophotometer (Thermo 
Scientific, Wilmington, DE, USA). A composite sample 
for pyrosequencing was prepared by pooling 200 ng of 
these purified PCR products of each sample. The pooled 
samples were submitted for pyrosequencing of the V3-V4 
region of the 16S rRNA gene on the 454 Life Sciences GS 
FLX platform using Titanium sequencing chemistry at 
GATC Biotech, Konstanz, Germany. 

16S rRNA gene sequence analysis 

Pyrosequencing data were analysed with a workflow 
based on QIIME vl.2, as described before [52]. Diver- 
sity metrics were calculated as implemented in 
QIIME vl.2. Hierarchical clustering of samples was 
performed using UPGMA with weighted UniFrac as a 
distance measure as implemented in QIIME vl.2. The 
Ribosomal Database Project classifier version 2.2 was 
performed for taxonomic classification [53]. The sig- 
nificance of the difference in relative abundance of 
specific taxa between males and females was calcu- 
lated using the Mann-Whitney U test as implemented 
in SciPy [54]. Additional data handling was done 
using in-house developed Python and Perl scripts. 

Statistical analysis 

Statistical analysis for the MA and microbiota data are 
described in the respective sections. Changes in DNA 
methylation were evaluated with an unpaired t test. For 
all tests, p values <0.01 were considered statistically 
significant. 

Results 

Gene expression in the small intestine and colon of 
prepubescent mice is sexually dimorphic 

To determine differences in gene expression in the intes- 
tine between male and female mice, whole -genome ex- 
pression profiling was performed on the small intestine 
(SI) and colon of 2-week-old male and female C57BL/6 
mice using microarray analysis (MA). Principal compo- 
nent analysis (PC A) was applied on the top 1,000 most 
variable genes present in the SI or the colon. The results 
obtained show distinct separation of the samples in two 
clusters by principal component (PC) 1, accounting for 
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90.6% of the expression variation. As seen in Figure 1A, 
one cluster contains the SI and the other cluster the 
colon samples. In addition, the samples within each 
cluster are divided by PC2 accounting for 3.3% of the 
variation, entirely separating the males and females in 
both intestinal segments. 

Next, genes showing significant (p < 0.01) sexually 
dimorphic gene expression were analysed in more de- 
tail for each of the two segments of the intestine sep- 
arately. Of the total number of 10,958 genes remaining 
after filtering, 275 genes displayed significant differen- 
tial expression in the SI including 98 male-dominant 
and 177 female-dominant genes (see Figure IB and 
Additional file 2). In the colon, a smaller number of 
genes (86) revealed sexually dimorphic expression pat- 
terns of which 58 genes demonstrated male-dominant 
and 28 genes showed female-dominant expression (see 
Figure IB and Additional file 3). Figure 1C shows that 
in total, seven genes displayed female-dominant ex- 
pression in the SI as well as in the colon. Five of them 
are located on the X chromosome and two are auto- 
somal genes. Of these seven genes, only Xist displays 
strong differential expression between the two sexes 
while for Kdm6a, Kdm5c, Pbdcl, Ddx3, Pisd-psi and 
Sptbn2 more subtle but significant differences were 
observed (Figure ID). The four genes displaying male- 
dominant expression in both segments of the intestine 
are KdmSd, Ddx3y, Uty and Eif2s3y. All of them are 
located on the Y chromosome and reveal, as expected, 
strong differential expression since these genes are not 
expressed by females. Furthermore, we identified two 
genes showing significant sexually dimorphic expres- 
sion in the two segments of the intestine but with an 
opposite sex effect. Cabp2 displays female-dominant 
expression in the SI and male-dominant expression in 
the colon while for Gm6581 male-dominant expression 
was found in the SI and female-dominant expression in 
the colon. In total, 349 different genes revealed sexually 
dimorphic expression in either the SI or the colon. 

By analysing the sexually dimorphic effects in the SI 
and colon in more detail, we found that the majority of 
the changes are relatively subtle in both segments of the 
intestine (see Figure IE, F). Evaluation of the fold changes 
(FC) of all sexually dimorphic genes in either the SI or the 
colon revealed that, in line with previously reported data 
in the liver, heart, adipose tissue and brain [19], the major- 
ity of the genes showed only small differences (FC <1.2) in 
gene expression between males and females (Figure 1G). 
The number of genes that differed more than 1.5 -fold be- 
tween males and females is very limited. 

To analyse the chromosomal distribution of the sexu- 
ally dimorphic genes, we calculated the chromosomal 
localization of the sexually dimorphic genes relative to 
the chromosomal localization of all 10,958 genes included 



in the analysis. As shown in Figure 1H, I, different distri- 
butions were found for the SI and colon, but in both seg- 
ments, the strongest effects were observed on the sex 
chromosomes (see Additional file 4). 

In summary, significant changes in gene expression 
were found in the SI and colon of prepubescent mice, 
with the highest FC between the sexes observed for 
genes located on the sex chromosomes. 

Inter-individual variation detected in many genes 
displaying sexually dimorphic expression 

Next, we concentrated our analysis on the genes display- 
ing the most pronounced sexually dimorphic effect in 
the SI and the colon. The top 25 genes displaying the 
strongest differential expression are presented in Table 1. 
The results reveal that, apart from the above-indicated 
four Y-chromosomal and X-chromosomal Xist genes 
showing strong sex-specific expression in both the SI 
and the colon, 20 genes are present, exhibiting a FC 
between 1.40 and 1.75. All these genes are located on 
autosomal chromosomes except for the X-chromosomal 
gene Kdm6a. Interestingly, for a number of these genes 
(indicated with §), intestinal expression and/or function- 
ing has been described [55-59] but until now sexually di- 
morphic expression has not been reported. 

By evaluating the expression levels of the top 25 
strongest sexually dimorphic genes in the individual 
mice in more detail, substantial variation in expression 
between the individual pups within the sex groups was 
observed. As seen in Figure 2A, mean expression values 
of Upklb, Dleu2, Fkbpll and Nts significantly differed 
between male and female mice, but at an individual 
level, some of the female mice show the same expres- 
sion level as a number of the male mice. The results for 
these genes are in contrast with the genes presented in 
Figure ID, where variation between the individual mice 
of one sex was extremely low. Hierarchical cluster ana- 
lysis (HCA) of expression levels of the subset of 329 
sexually dimorphic genes revealed separate profiles for 
male and female mice even in the absence of the genes 
located on the sex chromosomes. This result indicates 
that, despite the inter-individual variation in expression 
of the sexually dimorphic genes, male and female mouse 
pups reveal a clearly distinct expression profile (Figure 2B). 

Differential gene expression might have functional 
consequences for normal functioning and disease 
development between males and females 

Ingenuity pathway analysis (IPA) was applied to identify 
the biological functionality in physiology and disease as 
well as the functional networks of the sexually dimorphic 
genes. By analysing the physiological (Bio) functions in- 
depth, IPA revealed that the sexually dimorphic genes rep- 
resented various general cellular functions in both the SI 
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and the colon (Table 2). Regarding disease development, 
IPA revealed that in the SI, genes displaying sexually di- 
morphic expression were mainly involved in infection and 
inflammation, whereas in the colon also genes involved in 
cancer development were found to be differentially regu- 
lated (Table 2). 

In Table 3 the top networks having p scores of at least 
25 [60] are presented. The results obtained reveal that in 
the SI, sexually dimorphic expression is found for different 
subsets of genes involved in lipid metabolism (network 
identifier ID-1 and ID-2), genes playing a role in develop- 
mental networks (ID-3) and genes involved in "cell-to-cell 
signalling and interaction, haematological system develop- 
ment and function and immune cell trafficking" (ID-4). In 
the colon, two top networks with scores of at least 25 were 
found, one representing genes with a function in "car- 
diovascular system development and function, organ- 
ismal development and tissue morphology" (ID-1) and 
the other one containing genes with a function in 
"neurological disease, physiological disorders and pro- 
tein synthesis" (ID-2). 

Epigenetic mechanisms might contribute to differential 
expression of sexually dimorphic genes 

By analysing the total subset of 349 sexually dimorphic 
genes in more detail, we identified 9 genes exhibiting a 
function in epigenetic processes (see Table 4). Interest- 
ingly, apart from Mbdl, which is known for its binding 
to methylated DNA, all genes have histone-modifying 
activity. It can be speculated that the sexually dimorphic 
expression of these histone modifiers might contribute 
to the differential expression of a subset of the other 
sexually dimorphic genes. Since the amount of intestinal 
tissue of 2-week-old mouse pups is very limited and all 
available material was used for RNA and DNA isolation, 
we were not able to examine histone modifications in 
the SI and colon tissue in more detail. Since we could 
not measure histone modifications, we analysed DNA 
methylation as a surrogate marker of epigenetic regula- 
tion. DNA methylation has been shown to be involved 
in silencing of genes located on one of the two X chro- 
mosomes in somatic cells [61]. We measured DNA 
methylation by applying pyrosequencing on the promoter 
regions of a number of X-chromosomal genes. The results 
presented in Figure 3A, B clearly show significantly en- 
hanced methylation for all four CpGs analysed in the pro- 
moter of the Fundc2 promoter and for part of the five 
CpGs analysed in the Ddx3x promoter of female mouse 
pups. We used the top 25 of sexually dimorphic genes to 
evaluate the presence of promoter methylation as a poten- 
tial regulator of differential expression between the sexes. 
The promoter regions of the genes were identified by ap- 
plying Genomatix analysis, and the CpG content of each 
promoter was determined. As shown in Table 1, strong 



differences in the number of CpGs present in the differ- 
ent promoters were found, varying from a few CpG 
dinucleotides (i.e. Cyp2c66, Alox5ap, Cabp2, Mcptl, 
Mcpt2) to very CpG dense promoter regions (i.e. WT1, 
Syndigl, Dleu2, Kdrn6a). Pyrosequencing analysis re- 
vealed no differential methylation between males and 
females in any of the CpG sites studied that are present 
in the promoter region of Cyp2C66, AloxSap, Ccr3, 
Mcptl, Retnly and Nts, representing a subset of genes 
containing a relative low CpG content (see Figure 3C-H 
and Additional file 5). 

In summary, our data suggest potential involvement of 
mediators of histone modification in the regulation of 
sexually dimorphic effects, but we did not detect signifi- 
cant changes in CpG methylation in the promoter re- 
gions of a small subset of the top 25 sexually dimorphic 
genes between the two sexes. 

No sex effect on bacterial community compositions of the 
colon 

To evaluate the intestinal microbiota composition of 
the 2-week-old pups, deep sequencing of the 16S rRNA 
gene was applied on colonic luminal content collected 
from the mice. As seen in Figure 4A, the colon lumen 
of both male and female mice is predominantly colo- 
nized by Bacteroidetes and Firmicutes and, in addition, 
very small fractions of Actinobacteria and Proteobac- 
teria were measured. Evaluating the families colonizing 
the colon of the 2-week-old suckling mice revealed the 
presence of a variety of families in the colon lumen. 
The results presented in Figure 4B showed that hier- 
archical clustering using weighted UniFrac as a dis- 
tance measure separates the samples of all 15 mice 
included in the analysis into two major clusters, but 
it did not separate the males from the females. Re- 
dundancy analysis (RDA) revealed that the three lit- 
ter origins in which the mice were nursed during the 
2 weeks of their life have a dominant and highly sig- 
nificant effect on microbiota composition. This effect 
is visualized in the RDA plot presented in Figure 4C 
showing a clear separation in microbiota composition 
of the mice of the three litters whereby gender was 
included as a covariant. Multivariate statistics taking 
litter origin into account as a covariant showed a 
trend for sex effect on microbiota composition 
(Figure 4D), but this effect was not significant (p > 
0.01). Beta-diversity metrics calculated from the 
three different litters (Figure 4E) further confirmed 
the observation that litter but not sex (Figure 4F) has 
a strong effect on bacterial community compositions. 
Interestingly, for males, one specific operational taxo- 
nomic unit (OTU610 classified as Syntrophococcus) 
was identified as the most unique taxon associated 
with males (Figure 4G), which was confirmed by 
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(See figure on previous page.) 

Figure 1 Sexually dimorphic gene expression in the SI and colon of 2-week-old C57BL/6 mice. (A) PCA of the top 1,000 most variable 
genes present in the SI or the colon separates the SI from the colon by PC1 and males from females by PC2. (B) In the SI and the colon, 275 and 
86 genes displayed significant (p<0.01) sexually dimorphic expression, respectively. (C) Four of the significant sexually dimorphic genes revealed 
male-dominant expression in both the SI and colon while seven genes showed female-dominant expression in both intestinal segments. 
(D) Expression levels of the seven female-dominant genes in the SI and colon in all individual male and female mice. Correlation between male 
and female expression levels of the (E) 275 sexually dimorphic genes in the SI and (F) 86 sexually dimorphic genes in the colon. (G) The FC difference 
between males and females is very low for the majority of the sexually dimorphic genes. Chromosomal localization of the sexually dimorphic genes in 
the (H) SI and (I) colon relative to the total number of genes localized on each chromosome in the total selection of 10,958 included in the analysis. 
*p < 0.01, **p < 0.001, ***p < 0.0001 . 

\ J 



random forest analysis identifying it as the strongest 
predictor for males (see Additional file 6). Follow-up 
studies containing a larger sample size and isolated at 
different time points are required to validate and fur- 
ther explore this result. 

Taken together, our data reveal that litter origin has 
the strongest impact on bacterial community compos- 
ition and no significant effect of gender was observed. 



However, we cannot exclude that sex has an effect on 
the dominance of specific taxa. 

Discussion 

Sex-specific functionality and susceptibility to disease 
development of the intestine 

By applying MA analysis, we identified sexually dimorphic 
genes in both segments of the intestine of prepubescent 



Table 1 Top 25 genes displaying the strongest differential expression between males and females in the intestine 



Symbol 


SI 


Colon 


Chr 


Promoter length 


Number of CpGs 


Description 


Eif2s3y 


-49.68^ 


-46.96^ 


Y 


682 


21 


Eukaryotic translation initiation factor 2, subU 3, str. gene Y-linked 


Uty 


-48.21^ 


-44.48^ 


Y 


893 


28 


Ubiquitously transcribed tetratricopeptide repeat gene, Y chr 


Ddx3y 


-28.71 


-30.81 


Y 


643 


23 


DEAD (Asp-Glu-Ala-Asp) box polypeptide 3, Y-linked 


Kdm5d 


-15.67^ 


-21. 58^ 


Y 


757 


30 


Lysine (K)-specific demethylase 5D 


Upklb 


-1.12A 


-1.60^ 


16 


603 


11 


Uroplakin 1B 


Wtl 


-1.02 a 


-1.59^ 


2 


1282 


104 


Wilms tumour 1 homologue 


Syndigl 


-1.52<* 


-1.05A 


2 


603 + 1,225 


3 + 91 


Synapse differentiation inducing 1 


Nos2 


1.09* 


— 1 .41 ^ 


11 


673 


8 


Nitric oxide synthase 2, inducible 


Rab42 


-1.42^ 


-1.17A 


4 


783 


35 


RAB42, member RAS oncogene family 


Nts § 


—1 .41 ^ 


-1.15A 


10 


794 


6 


Neurotensin 


Cyp2c66 


1.45? 


-1.07A 


19 


613 


1 


Cytochrome P450, family 2, subfamily c, polypeptide 66 


Mmp25 


1.43? 


1.05 a 


17 


893 


36 


Matrix metallopeptidase 25 


Slc7a11 § 


1.40? 


1.05 a 


3 


829 


18 


Solute carrier family 7, member 1 1 


Nucb2 § 


1.46? 


1.07 a 


7 


666 


36 


Nucleobindin 2 


Alox5ap § 


1.44? 


1.13A 


5 


705 


4 


Arachidonate 5-lipoxygenase-activating protein 


Gm 12696 


-1.13A 


1.41? 


4 






Predicted gene 12696 


Cabp2 


1.50? 


-1.06 a 


19 


603 + 820 


3 + 8 


Calcium-binding protein 2 


Fkbpl 1 


1.56? 


1.09A 


15 


1041 


27 


FK506-binding protein 1 1 


Mcptl 


-1.02 a 


1.54? 


14 


606 


5 


Mast cell protease 1 


Mcpt2 


-1.12A 


1.66? 


14 


601 


1 


Mast cell protease 2 


Dleu2 


— 1.1 1 A 


1.71? 


14 


798 + 601 


84 + 45 


Deleted in lymphocytic leukaemia, 2 


Ccr3 


1.59? 


1.50A 


9 


730 


3 


Chemokine (C-C motif) receptor 3 


Retnly § 


1.72? 


1.52A 


16 


638 


4 


Resistin-like gamma 


Kdm6a 


1.59? 


1.57? 


X 


693 + 601 


51 +69 


Lysine (K)-specific demethylase 6A 


Xist 


104.49? 


182.91? 


X 


1,065 


29 


Inactive X-specific transcripts 



^Significant (^-dominant expression FC >1.4. 
^Significant $ -dominant expression FC >1.4. 
A FC without a significant sexually dimorphic effect. 

§ Previous studies revealed an intestine-specific function and/or expression of this gene. 
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Figure 2 Inter-individual variation in gene expression within the groups of male and female mice. (A) Expression levels of Upkbl, Dleu2, 
Fkbpll and Nts in the SI and colon of the individual male and female mice. (B) Hierarchical clustering of the expression levels of the subset of 
329 sexually dimorphic genes revealed separation of the males and females in two distinct clusters. *p < 0.01. 



C57BL/6 mice, with the most abundant effect observed in 
the SI. The vast majority of genes displayed only small 
expression differences between males and females. This 
effect is in line with previously reported sexually di- 
morphic effects found in adult liver, brain, muscle and 
adipose tissue of 334 mice [19]. A small subset of genes 
displayed more pronounced sexually dimorphic effects 
including Nts, Nucb2, Slc7all, AloxSap and Retnly, 
which are genes that have previously been linked to 
intestinal expression and/or functioning. Our results 
reveal male-dominant expression of Nts, a gene expressed 
by the enteroendocrine N cells in the intestine. The 13- 



amino acid gene product is excreted and serves multiple 
functions including stimulating pancreas and biliary secre- 
tion [62,63] and affecting gastric, SI and colonic motility 
[59,64]. Nucb2, which displays female-dominant expres- 
sion in the SI, is expressed at different positions of the 
gastrointestinal system including the duodenum [65]. 
Nefstatin-1 is the 82-amino acid excreted peptide derived 
from Nucb2 and has been shown to play a role in food in- 
take regulation by inducing an anorexic signal in the 
hypothalamus [66]. Slc7all mediates the transfer of 
amino acids over the cellular membrane, and its ex- 
pression in the intestine has been described before [56]. 



Table 2 Top Bio functions and diseases linked to the sexually dimorphic genes 

Bio functions p value #mol Diseases p value #mol 

SI 



Cellular movement 


8.04E 


-14 to 6.18E 


-3 


67 


Inflammatory response 


1.23E 


-12 to 6.1E 


-3 


60 


Cellular function and maintenance 


5.95E 


-11 to 6.1 E- 


-3 


68 


Hypersensitivity response 


7.04E 


-8 to 9.51 E 


-4 


24 


Cellular compromise 


3.84E 


- 9 to 2.44E - 


-3 


23 


Infectious disease 


6.1 1E 


-5 to 5.59E 


-3 


29 


Cell signalling 


2.01 E 


- 8 to 5.71 E- 


-3 


40 


Hereditary disorder 


1.16E 


-4 to 2.78E 


-3 


6 


Molecular transport 


2.01 E 


-8 to 5.71 E- 


-3 


62 


Inflammatory disease 


1.16E 


-4 to 3.29E 


-3 


43 


Colon 




















Cellular development 


1.12E 


-4 to 4.19E- 


-2 


19 


Haematological disease 


6.06E 


-4 to 3.36E 


-2 


11 


Cellular growth and proliferation 


1.12E 


-4 to 4.19E- 


-2 


23 


Inflammatory response 


1.5E- 


3 to 4.19E- 


2 


12 


Cell morphology 


1.4E- 


4 to 4.31 E- 


2 


17 


Cancer 


1.57E 


-3 to 4.21 E 


-2 


19 



#mol number of molecules included in the indicated Bio function or disease. 
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Table 3 Top networks of sexually dimorphic genes 



ID Top networks mol/score a 

si S vs ? 

1 Lipid metabolism, small-molecule biochemistry, vitamin and mineral metabolism 24/41 
ABCG1$,ALDH1A3$,ALS2CR12cT,CD300 

mir-1 81 $,MYBL1 c?,MYBL2$,NCAM2cJ,PREX1 $,S1 00A1 cT,SELPLG$,SH3BP1 $,SLC46A3cJ,TIFAcJ,TOLLIPcJ,UBA3cJ 

2 Lipid metabolism, small-molecule biochemistry, haematological system development and function 18/28 
ACER3cTASAH1c?,BTK$,CD22$,LAPTM5$,NFATC1$,P2RX7$,PIK3AP1$,PIK3CD$,PLA2G7$,PLA2G7$,PLA2G10$, 
PLCG2$,PLD2$,PRR7$,RASGRP4$,TNRcJ,VPS45cJ 

3 Haematological system development and function, tissue development, cellular movement 17/26 
AHR$,CCL1 1 $,GNPDA1 c?,HCLS1 $,IL1 5c?,KDM5Dc?,KDM6A$,LAMP1 c?,MBD1 cT,MZB1 $,P2RX1 $,PLCB2$,PYGL$, 
RAC2$,SRGN$,TRIM5c?,UtycJ 

4 Cell-to-cell signalling and interaction, haematological system development and function, immune cell trafficking 17/26 
BET1cT,BMP7cT,Ceacam1/Ceacam2^ 

RUNX3$,SAMHD1 $,SH3KBP1 $,SH3PXD2BcJ,STAP1 $,TRAM2c? 

Colon S vs $ 

1 Cardiovascular system development and function, organismal development, tissue morphology 20/44 
CABP4cJ,DDX3X$,Dleu2$,DLK1 c?JTGA2Bc?,KLK3$,MAP3K8c?,mir-1 0c?,mir-1 81 c?,MTCH2$,NOS2cJ,NPR3cJ, 
PDGMCcT,PRKCZcT,PRXcT,STX6cT,SYT9cTjNMRSM19cTTNMRSM11^ 

2 Neurological disease, psychological disorders, protein synthesis 16/33 
BIVm$,CWC22$,EIF2S3 c T / FUNDC2$ / KDM^ 

STAG3cT,SURF6$,Utyc?,WDR44$ 
Networks with p scores >25 are presented. 

5 genes displaying ^-dominant expression, $ genes displaying $ -dominant expression. 
a Focus molecules/score. 



Our study revealed female-dominant expression of 
Slc7all in the SI of the prepubescent mice. AloxSap is 
involved in the synthesis of leukotrienes which are 
lipid-signalling molecules derived from arachidonic 
acid and showed higher expression in the SI of female 
mice. Leukotrienes are known to initiate and amplify 
inflammation, and expression of AloxSap in intestinal 
cells has previously been reported [57]. Finally, Retnly 
has been reported to be expressed in the intestine 



and has been suggested to play a role in intestinal 
tract-mediated regulation of insulin sensitivity [67] 
and in inflammatory processes [68]. Retnly showed 
significant female-dominant expression in the SI of 
the 2-week-old mice. In addition, the list of genes 
revealing more subtle changes in gene expression be- 
tween males and females included additional trans- 
porters and genes involved in lipid metabolism. 
In conclusion, our results indicate that sexually 



Table 4 Genes displaying sexually dimorphic expression and exhibiting a role in epigenetic processes 



Symbol 


FC-SI 


FC-Colon 


Description 


Chr 


Epigenetic activity 


Mbd1 


-1.13^ 


-1.03 


Methyl-CpG-binding domain protein 1 


18 


Binding to methylated DNA 


Phf20 


1.02 


-1.14^ 


PHD finger protein 20 


2 


Histone acetylation 


Uty 


-48.21^ 


-44.48^ 


Ubiqu. transcribed tetratricopeptide repeat gene, Y chr 


Y 


Histone demethylation 


Kdm5d 


-15.67^ 


-21.58^ 


Lysine (K)-specific demethylase 5D 


Y 


Histone demethylation 


Kdm6a 


1.59? 


1.57? 


Lysine (K)-specific demethylase 6A 


X 


Histone demethylation 


Kdm5c 


1.32? 


1.28? 


Lysine (K)-specific demethylase 5C 


X 


Histone demethylation 


Bcorll 


1.15? 


1.04 


BCL6 co-repressor-like 1 


X 


Associates with HDAC activity 


Hdacl 1 


1.14? 


-1.11 


Histone deacetylase 1 1 


6 


Histone deacetylation 


Piwil4 


1.17? 


-1.03 


Piwi-like RNA-mediated gene silencing 4 


9 


Histone methylation 



Significant (^-dominant expression; ^significant ^-dominant expression. 
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Figure 3 DNA methylation analysis of the promoter regions of genes displaying sexually dimorphic expression. Significant differential 
methylation was found in the promoter region of two X-chromosomal genes: (A) Funcd2 and (B) Ddx3x. No differential CpG methylation was found in the 
promoter region of the autosomal genes (C) Cyp2c66, (D) Alox5ap, (E) Ccr3, (F) Mcptl, (G) Retnly and (H) Nts, all displaying significant sexually dimorphic 
expression with a FC difference >1 .4 between males and females. Blue bars: males; red bars: females. *p < 0.01, **p < 0.001, ***p < 0.0001. 



dimorphic expression of genes might induce distinct 
functioning of both segments of the intestine in males 
and females. 

The results obtained by IPA carried out on the signifi- 
cantly differentially expressed genes imply that also dis- 
ease development might be regulated differently between 
the two sexes. Predominantly, genes involved in infec- 
tious disease and inflammatory processes were differen- 
tially expressed between the two sexes and— again— this 
effect was most pronounced in the SI. In addition, the 



top 25 list of genes exhibiting the most pronounced 
sexually dimorphic expression in either the SI or colon 
included four genes (Nos2, Mcptl, Mcpt2, Ccr3) that 
carry out functions in immune response and inflamma- 
tion. In humans, one of the most frequently occurring 
diseases of the intestine is inflammatory bowel disease 
(IBD), which includes two chronic diseases that cause 
intestinal inflammation: Crohn's disease and ulcerative 
colitis (UC). Both are multifactorial disorders and factors 
playing a causal role in disease development are genetic 
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Figure 4 (See legend on next page.) 
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(See figure on previous page.) 

Figure 4 Microbiota composition analysis of the colonic luminal content revealed no significant differences between male and female 
prepubescent mice. (A) Bacteriodetes and Firmicutes dominate the colon lumen in 2-week-old male and female pups. (B) Hierarchical clustering 
using weighted UniFrac as a distance measure of the male and female mouse pups revealed no clustering of the sexes. (C) RDA plots show 
microbiota differences based on litter origin, but (D) only marginal and not significant differences were found between male and female mice. 
Boxplot diagrams of the beta-diversity showed (E) a more distinct difference of the mice of litter 1 compared to litters 2 and 3 but (F) no 
significant differences between males and females. (G) OTU610 (classified as Syntrophococcus) was identified to being most uniquely associated 
and hence predictive for the male sex. 



predisposition, the intestinal microbiota, the immune 
system and environmental factors— in particular nutrition 
[69-71]. There is currently no consensus regarding the ef- 
fect of sex on disease development of both Crohn's disease 
and UC. While some of the epidemiological studies did 
observe sexually dimorphic effects, other studies did not 
[72]. Interestingly however, several of the identified sexu- 
ally dimorphic genes containing a function in inflamma- 
tory processes have previously been linked to IBD 
including Ccr3 [73], Cclll [74], Ncfl [75] and Tnfr [76]. 
Apart from genes playing a role in immune response and 
inflammation, we identified cancer-relevant genes in the 
colon showing sex-differential expression using IPA. Fur- 
thermore, the top 25 of most differentially expressed genes 
includes two genes that have previously been linked to 
colorectal cancer: Mmp25 [77] and WT1 [78]. Taken to- 
gether, our data reveal that genes involved in disease de- 
velopment of the intestine in adults show differential 
expression between the sexes. It should be taken into ac- 
count that in this study, we analysed prepubescent mice. 
Samples collected at other ages, as well as human samples, 
are required to further explore the molecular mechanisms 
that might be responsible for differences in intestinal dis- 
ease development between males and females. 

Inter-individual variation within the sexes 

Even though the results of the MA analysis revealed signifi- 
cant sex differences for the expression of numerous genes, 
most fold changes are relatively low. It might be anticipated 
that, when a larger sample number would have been ana- 
lysed, the number of genes revealing significant sexually di- 
morphic expression might have been higher. Contributing 
to the relatively low fold changes, we found substantial 
inter-individual variation between the individual mice 
within one sex, especially for a subset of genes displaying 
stronger expression differences between males and females 
(FC 1.4-1.75). Although C57BL/6 is an inbred mouse 
strain, strong phenotypic variations between male mice at 
different ages in response to exposure to a high-fat [79] or 
medium-fat (Rusli et al: A weekly alternating diet between 
caloric restriction and medium-fat protects the liver from 
development of NAFLD in middle-aged C57BL/6J mice, 
submitted) diet have been described previously by us. 
Nonetheless, hierarchical clustering of the significantly dif- 
ferentially expressed genes— after exclusion of genes located 



on the sex chromosomes— still separated the male and fe- 
male mice in two distinct groups. This result indicates that 
male and female expression profiles of the sexually di- 
morphic genes are distinct. This effect is independent of 
the variation in expression within the two sexes and is de- 
tectable even in the absence of sex-linked genes. Studies in- 
cluding larger numbers of mice are required to evaluate 
inter-individual variation of the sexually dimorphic genes in 
more detail. 

Sexually dimorphic effects pre- and post-puberty 

In our study, we analysed sexually dimorphic effects in 
the intestine of prepubescent mice and found that in 
these 2-week-old mice pups, already substantial changes 
in gene expression were detectable in particular in the 
SI. We have previously reported physiological and mo- 
lecular differences in the liver between the two sexes of 
these young mice in response to maternal exposure to a 
Western-style diet [43]. By comparing the significantly 
sexually dimorphic genes between the intestine and the 
liver, we found that most of the genes displaying differ- 
ential expression between the sexes in the two segments 
of the intestine were also found in the liver. However, 
overall the overlap in genes displaying significant sexu- 
ally dimorphic expression between the liver and the in- 
testine was very low (Additional file 7). By comparing 
differential expression of male and female adult mice in 
different tissues, Yang and colleagues observed high tis- 
sue specificity for the sexually dimorphic genes in the 
liver, brain, adipose tissue and muscle [19]. By analysing 
the common functional categories of the sexually di- 
morphic genes, they found that only steroid and lipid 
metabolism were shared between the different tissues. 
Interestingly, we also found differential expression of 
genes contributing to lipid networks in the SI. The rea- 
son why we did not observe an effect in steroid metabol- 
ism might be due to the fact that this study has been 
carried out in 2-week-old pups when hormonal differ- 
ences between the sexes are still extremely low. Evalu- 
ation of the 349 sexually dimorphic genes revealed that 
this list did not include well-established oestrogen- or 
androgen-responsive genes. 

Our analysis included samples collected of prepubes- 
cent mice. In previous studies of Conforto and Waxman 
[26] and Kwekel et al. [21], sexually dimorphic gene 
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expression was analysed in pre- as well as post-pubertal 
animals. Both studies revealed that, although sexually di- 
morphic effects are detectable in pre-pubertal animals, 
the effects are much more pronounced in the adult 
phase of life. These results imply that also the difference 
between males and females for the two segments of the 
intestine might be much stronger in adult mice. On the 
other hand, Kwekel and coworkers noticed that part of 
the sexually dimorphic effect detected in the youngest 
animals disappeared at older ages [21], indicating differ- 
ences in the kinetics of the developmental processes. 
Interestingly, according to data published by Mank and 
colleagues, the kinetics over the life cycle of sexually di- 
morphic genes in chicken differ between males and fe- 
males [24]. Whether similar effects occur in the intestine 
needs to be explored in further studies. 

Variation in microbiota composition between the mice is 
driven by litter origin but not by sex 

The newborn gut is essentially sterile. Inoculation and 
colonization of the colon starts already during delivery, 
and time is required to allow the first settlers to establish 
a dense microbial population. The final population col- 
onizing the colon includes thousands of bacterial species 
that belong to a limited number of phyla. In the adult in- 
testinal lumen, Bacteroidetes and Firmicutes are the 
dominant phyla, but during the colonization phase, 
bifidobacteria species (members of the Actinobacteria 
phylum) are usually dominant in breastfed infants [37]. 
However, in our 2-week-old mouse pups, an extremely 
limited amount of bifidobacteria was observed and the 
colon microbiota was almost exclusively made out of 
Bacteroidetes and Firmicutes. Our results indicated that 
variation between bacterial communities was dominated 
by litter origin while no significant difference between 
males and females was found. Most likely, the micro- 
biome of their mothers is a critical component of the 
offsprings' microbiome at this early pre-weaning phase 
and might dominate over other potential modifiers. 
However, we cannot exclude that sex has an effect on 
the dominance of specific taxa. These observations are 
consistent with previous reports, which showed that gut 
microbiota differences become evident in adult male and 
female mice and are driven by androgen hormones, but 
are absent in prepubescent mice [80,81]. Furthermore, 
various other influencing factors including cage and ma- 
ternal effects have been reported previously [82,83]. 

Localization on the sex chromosomes together with 
epigenetic effects might be responsible for the sexually 
dimorphic gene expression 

As indicated above, sexually dimorphic gene expression 
is much more robust in the adult phase of life than 
during the prepubescent phase due to circulating sex 



hormones. Although very low levels of sex hormones 
are already present in early life, there is increasing un- 
derstanding that also hormone-independent pathways 
of sexual differentiation exist [84]. Alternative mecha- 
nisms responsible for molecular regulation of sexually 
dimorphic effects are (1) expression of genes located 
on the sex chromosomes and (2) epigenetic differences 
between the sexes. Our study revealed relative abun- 
dance of sex chromosomal localization amongst the 
sexually dimorphic genes. Of the total subset of 349 
intestinal sexually dimorphic genes, 19 were located on 
the X chromosome and 4 on the Y chromosome. Fur- 
thermore, our results showed that, consistent with 
previously reported data [85], the majority of the X- 
chromosomal genes displayed female-dominant expres- 
sion and apart from the well-known escapees Xist, Kdm6a 
and KdmSc, only marginal expression differences between 
males and females were detected. The Y-linked copy of 
Kdm6a and Kdm5c, Uty and Kdm5d, respectively [8], re- 
vealed strong expression in the intestine of the 2-week- 
old mice pups. Many of the sexually dimorphic genes 
located on the sex chromosomes display a regulating 
function in transcription (i.e. Btk, Ddx3x, BcorLl, Uty) 
or translation (Eif2s3y), and it can be speculated that 
these genes contribute to the differential expression of 
the autosomal genes. 

Other potential mechanisms responsible for sexually 
dimorphic expression are epigenetic modifications. Pre- 
vious studies have provided evidence that changes in his- 
tone modifications [22,42] as well as alterations in DNA 
methylation [23] might contribute to differential gene 
expression between the two sexes. Our study displayed 
differential expression of a subset of genes involved in 
regulating histone modifications supporting the previ- 
ously reported results. By analysing the promoter regions 
of the top 25 strongest sexually dimorphic genes, we ob- 
served high diversity in the CpG content ranging from 
extremely high CpG dense to just one or a few CpGs. 
We analysed DNA methylation of various genes contain- 
ing low or medium CpG dense promoter regions and 
found no significant changes in DNA methylation in any 
of the 14 analysed CpG sites between males and females. 
Weber et al. identified a link of CpG density of promoter 
regions with their methylation status [86] which might 
be responsible for our findings instead of sex- differential 
effects. Moreover, pyrosequencing has the limitation to 
analyse only selected CpGs/regions and is not a genome- 
wide approach. Therefore, differential DNA methylation 
might occur on CpGs in the promoter regions that have 
not been included in our analysis. Furthermore, changes 
in gene expression might also be regulated by changing 
the methylation status of CpGs present in enhancer re- 
gions or in the gene body of the genes, which were not in- 
cluded in our analysis. Genome-wide DNA methylation 
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analysis to investigate CpG methylation on a more exten- 
sive scale is required to determine the involvement of 
DNA methylation as underlying regulatory mechanisms of 
sexually dimorphic gene expression in the SI and colon. 

Conclusions 

Our study revealed sexually dimorphic genes in both 
segments of the intestine. Functional analyses of these 
genes pointed towards differences in normal physio- 
logical functioning as well as disease development be- 
tween males and females. No significant differences in 
bacterial community composition were found, although 
sex-specific taxa might exist. Since prepubescent animals 
have been used in this study, we can conclude that, even 
in the absence of circulating sex hormones and in the 
presence of identical microbiota in the colon lumen, 
there is intrinsic sex-specific gene regulation. Since the 
intestine fulfils a central role in whole-body health, the 
observed molecular sexually dimorphic effects might 
contribute to the differences in basic physiology, body 
composition and susceptibility to and progression of a 
broad variety of non-communicable diseases between 
males and females. 
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